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DETEEMINATION OF ROTOR HARMONIC BLADE 
LOADS FROM ACOUSTIC MEASUREMENTS* 


1,0 INTRODUCTION 

The sovmd radiated by a rotating blade is characterized by a combination of 
broadband noise and discrete frequency noise. The noise spectra generated 
by the main rotor of a helicopter during normal operation are generally domi- 
nated by discrete frequencies below 500 Hz, 


The radiation of discrete tones or rotational noise results from the cyclic 
interaction of aerodynamic forces on the acoustic media. The theoretical 
derivation of rotational noise due to a steady uniform force over a rotor disk 
is contained in a classic report by Gutin (Reference 1), In recent years, it 
has been recognized that magnitudes of the sound harmonics are primarily 
controlled by the presence of a nonuniform distribution of aerodynamic forces 
over the rotor disk. When expressed in terms of a Fourier series in angular 
position about the rotor axis, the nonuniform forces are commonly referred 
to as blade-load harmonics. 


The magnitudes and frequency characteristics of the rotor harmonic blade 
loads have been studied experimentally by measuring the fluctuating pressures 
on the rotor blades with surface mounted pressure transducers. A recent 
report by Hosier and Ramakrishnan (Reference 2) contains a discussion of 
the experimental procedures and examples of typical surface pressure data. 
This direct measurement approach poses two major difficulties. One is the 
inherent experimental problem of recovering an electrical signal from a 
transducer rotating with the blade. The other difficulty relates to the con- 
ceptual problem of interpreting the net effect of measurements taken at iso- 
lated points along the blade. As discussed in Reference 2, the accuracy of 
rotational noise level predictions are dependent on the interpreted blade 
loading distribution. 


The approach taken in this investigation was to deduce the blade- load har- 
monic amplitudes from measurements of the discrete frequency acoustic 
levels taken at a position on the rotor axis. The advantages of this approach 


*The contract research effort which has led to the results in this report was 
jointly sponsored by NASA-Langley Research Center and USAAMRDL 
(Langley Directorate), 



were that it required a less complicated measurement procedure and offers a 
direct means of determining the net effective harmonic loading. 


The initial phase of the program was to undertake a review of current rota- 
tional noise theories. Theoretical approaches applicable to the blade meth- 
odology, presented in References 3 through 6 were reviewed in detail. The 
approach of Lows on and Ollerhead (Reference 3) consisted of modeling the 
acoustic properties of the rotor blade by means of a rotating point force. A 
solution to the three-dimensional wave equation was developed, expressing 
the radiated sound in terms of the harmonic blade-load components. The 
approaches of both Wright (References 4 and 5) and Morse and Ingard (Ref- 
erence 6) consisted of modeling the forces applied to the acoustic media in 
terms of a distributed loading over the rotor disk. Wright utilized a formu- 
lation as developed by Lighthill to express the sound radiated by a dipole 
force distribution, Morse and Ingard took the mathematical approach of using 
a Green's function to relate the radiated sound pressures to blade-load forces. 
The solutions of both Wright, and Morse and Ingard can be compared with 
that of Lows on and Ollerhead by considering the source distribution as a con- 
centrated loading at a given effective rotor radius. A comparison of the 
analytic expressions derived from these three approaches showed basic 
functional agreement. 


Section 2. 0 of this report provides a general presentation of the rotational 
noise theory utilized for calculations and predictions contained in the report. 
The theoretical section is divided into two parts. The first part provides the 
derivation of the relationship between on-axis acoustic harmonics and cor- 
responding blade-load harmonics. The second part contains a derivation of 
the general rotational-noise prediction model using an approach similar to 
that of Reference 6. 


Section 3. 0 of the report provides a compilation of the blade-load harmonic 
values derived from acoustic measurements taken at the NASA-Langley 
Research Center for the following configurations: 

a. OH-6 helicopter in hover at 24.4 meters (80 feet). 

b. Bell 205 helicopter in hover at 24, 4 meters (80 feet). 

c. NASA LRC helicopter rotor test tower with a type UHl-D rotor blade. 

Section 3.0 also contains a comparison of the measured off-axis sound har- 
monics with predictions based on the derived blade-load harmonic values. 

The conclusions of this investigation are contained in Section 4. 0. 
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2, 0 ROTATIONAL NOISE THEORY 


As discussed in References 1 through 6, the acoustic source mechanism 
created by the interaction of a rotor blade with a fluid medium is dipole in 
nature. The aerodynamic forces within the rotor disk are described in 
terms of a continuous distribution of acoustic dipoles. The sound field due 
to this source distribution is derived by using the mathematical properties of 
a Green's function. As described in Reference 6, the Green's function cor- 
responds to an expression for the acoustic pressure at an observation position 
due to a point source of unit strength located at a given source position. The 
total sound pressure due to a distributed acoustic source can then be deter- 
mined by integrating the product of the Green's function and a source dis- 
tribution function. To maintain a degree of simplicity in the mathematical 
model, the aerodynamic force distribution along the span of the rotpr is 
_taken to act at a single point. The model does, however, permit a consider- 
ation of the effect of loading variations across the chord of the rotor blade. 


The theoretical treatment is applied first to derive the relationship between 
the near-field on-axis sound levels and the rotor blade-load thrust harmonics. 
Secondly, the general relationship between far-field sound harmonics and 
blade-load harmonics is derived. 


2 , 1 Calculation of Harmonic Blade Loads from Axial Sound Harmonics 

The radiation of discrete frequency sound from a rotating blade is a physical 
phenomenon exhibiting a complex dependence on phase. The following theo- 
retical development models the effect of the rotating blades as a distribution 
of dipole forces acting on the plane surface enclosed by the blade tips. The 
sound pressure at an arbitrary observation point is the instantaneous sum of 
contributions from all points of the force distribution, each with appropriate 
phase. Mathematically, the expression for sound radiated along the rotor 
axis is simplified from the general expression for radiated sound. This 
simplification results because all sound radiated from points at a fixed blade 
radius arrive in phase at any observation point on the rotor axis. The forces 
acting on the rotor disk surface are taken to act at a given effective radius. 


A diagram of the cylindrical coordinate system to be used for the derivation 
of the harmonic blade loads from axial sound harmonics is shown in Figure 1. 
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The forces acting on the rotor disk surface are taken to act at a single 
effective radius. The positive z-axis is the direction of blade thrust, and the 
blade is rotating in the direction of positive uu. Besides the approximations 
inherent in the point force concept, the analysis is otherwise exact in the 
classical sense, and therefore, provides both near-field and far-field 
results. 


As discussed in Reference 7, the Green's function mathematically represents 
a solution to the Helmholtz equation expressing the spatial factor of a wave 
generated by a unit simple -harmonic point source. Using the coordinates 


(Eq. 1) 


of Figure 1 , 

the Green 

's function for sound radiated along the rotor 

given by 





exp (iskR) 


g 

41TR 


2 

2 2 

where 

R 

z + r 

o 


z = 

height of observation point above rotor 


r = 

radial distance 


= wave number of fundamental sound harmonic 


= harmonic integer 


The Green's function for sound radiated by a point dipole source is derived by 
taking the gradient of the simple point source function. For the force com- 
ponent in the z direction, the dipole Green's function thus becomes 


3 

z 9 z 


iskz iskR z iskR 
e - e 

4TTR 4TTR 


(Eq. 2) 


Since the Green's function for sound along the axis is independent of 0, the 
contribution to sound radiation due to the drag component in the plane of the 
rotor disk is zero. That is 


g 


0 


1 


r 

o 



(Eq. 3) 
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The next item required is an expression for the distribution of force in the 

z direction over the rotor disk. As discussed in references 3 through 6, the 

force distribution over the rotor disk may be represented by a double Fourier 

series expressing periodic conditions in both time and in azimuthal angle 

0 , Using complex notation, the Fourier series can be put in the form 
o 

CO 

^o “ ^(^o ~ (Eq. 4) 

s,q=~ 


where 


thrust distribution 


II 

average thrust force 
length of circumference 


a = 

s 

normalized chord load coefficient 


Pq = 

normalized blade thrust load harmonic coefficient 


s = 

sound harmonic number 


B = number of blades 

q = blade load harmonic number 


6( r —a) 
o 


Dirac delta function 


a = effective radius 


In general, the coefficients ttg and Pq are complex numbers. The coefficient 
CLs represents the normalized harmonic force component at a fixed point on 
the rotor disk due to the periodic passage of a rotor blade. The value of the 
coefficient ttg for a given harmonic number is a function of the rotor blade 
chord width, the blade repetition rate, and the chord loading distribution. 

At low soxand harmonic frequencies where the acoustic wavelengths are much 
larger than the chord width, the coefficients approach unity. For the range 
of blade-load harmonics considered in this study, the effect of chord loading 
is neglected, and all the coefficients represented by a are set equal to unity. 
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The blade -loading coefficient 6q represents the periodic force component 
that is experienced at a position on the blade as it rotates through a 
time -independent, but nonuniform flow. The normalization of 8q is such that 
the steady or net thrust load is equal to unity, and that each Pq with o 
represents the ratio of the qth harmonic load to the steady load. 


Using the mathematical properties of the Green's function, the acoustic 
pressure along the rotor axis is determined by integrating the product of the 
Green's function and the source distribution over the rotor disk, with radius 


P = 


2TT r 



T 


f„e r dr d0 
T®z o o ^o 


(Eq. 5) 


Replacing f-j- and g^ with the previously indicated expressions gives the 
equation 


P = F. 


^exp (- 


s, q= 


3uut ) a B 

' s q 


2TT r 

/ / 

r\ r\ 


&(i 


o 


iskz z 


i ( sB+q ) 0 + iskR 

exp 

2 3 

Uttr 4TTR j 




r, dr d0 
o o o 


(Eq. 6) 


Next, using the property of the Dirac delta function, the integration over the 
radius of the rotor disk is given by 


P = aP 


X a B 

s q 


00 

^ ^ j^exp|iskR^— isuut) 


s, q=-“ 
iskz 


2 3 

4TTR 4rTR J 
a a 


2rr 

j exp ( i ( sB+q ) 0^ j d0. 


(Eq. 7) 


with R 

a 



2 

+ a 
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Considerable simplification results from the integration. All terms in 
the double series are equal to zero except those in which sB equals -q. The 
remaining terms increase in amplitude by a factor of ZTT: 


P = 2lTaF° 
T 



iskz 

-4TrR^ 

a 


z 

4TTR^ 

a 




(Eq. 8) 


To examine the physical significance of this expression, it is helpful to con- 
vert from complex to real notation. This is readily accomplished by first 
expressing PsB ^ terms of an amplitude and a phase. 


That is 


®sB ‘ 


By convention, the value p is similarly expressed as 

- s B 


(Eq. 9a) 


®-sB = 


(Eq. 9b) 


If by some circumstance PsB were a real quantity, it would represent a har- 
monic blade load which was symmetric about some diameter in the rotor 
disk. Conversely, if ■^gg were purely imaginary, it would represent an 
antisymmetric blade load. To close approximation, the coefficients Ctg 
represent a symmetric time dependence and will be taken as real values. 

The summation limits of Equation 8 are changed using the relationship 





(Eq. 10) 


Rearranging terms and converting from exponential to trigonometric forms 
provides 


P = ZTTaF^ ^ ^ 2a^ 1 | | ^ ^ cos v| 


s= 1 


4TTR 


4TTR' 


where v = skR — suut — Y ^ 
a sB 


(Eq. 11) 
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The s = 0 term in the 'expression for P has been dropped since it corresponds 
to a DC pressure wave and does not contribute to the acoustic pressure. The 
coefficients of the sine and cosine terms in this series represent the far -field 
and near-field terms, respectively. 


In practice, the quantities ultimately measured are the harmonic values of 
the mean-square pressure. Taking the mean-square value of the above series, 
equation 1 1 yields the following 


P 


2 




1 1 2 

f/ skz \2 

/ z \ 2 I 

lPsB> 

2 + 3 

L 4TTR ' 

'4TTR ' -J 


a a 


(Eq. 12) 


where 


total mean 


thrust = 2TTaF^ 


prms 

s 


root-mean- square magnitude of acoustic harmonic 
component of harmonic number s 


By equating the above series term-by-term, the magnitudes of the blade-load 
harmonics are solved in terms of the individual acoustic harmonics. 

That is 


rms 


sB 


/TVs 


( skz 
41TR^ 


2 



1 - 1/2 


(Eq. 13) 


For the special case of sound radiation along the rotor axis, each acoustic 
harmonic component is related to an individual blade-load harmonic. It is 
important to note from Equation 13 that the acoustic harmonic designated 
by integer s is related solely to the blade harmonic with harmonic number 
sB. Taking the example of a rotor with two blades, on-axis measurements 
of the acoustic harmonics could be used to determine only the even-numbered 
blade- loading harmonics. Also of importance, in the above expression, is the 
result that the relationship between on-axis acoustic harmonic and blade-load 
harmonics is independent of the blade-load harmonic phases. 
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2. 2 D erivation of General Analytic Model 

The theoretical development for the general case of sound radiation from a 
rotor disk is an extension of the analytic technique used in the previous 
section. In view of the greater complexity, it is necessary to make pertinent 
approximations to arrive at a tractable solution. The coordinate geometry 
used for the derivation is shown in Figure 2. 


Proceeding as before, the initial step is to develop a functional form for the 
Green's function. As detailed in Reference 6, the contribution at any 
observation point due to a unit point source in the plane of the rotor disk is 
given by 


g = 


41T 



(Eq. 14) 


where the vectors'? and represent the spherical coordinate points (r, 9, 0) 
and (tq, 0o)» respectively. By limiting our considerations to the far- 
field case where r » Tq, the distancel"? - ?"o| is approximated by the 
following expression 


r - 


r ~ r — r sin 
o — o 


cos 


( 0-0 ) 


(Eq. 15) 


Accordingly, the Green's function is then approximated by 


g = 


1 

4TTr 


exp 



iskr 

o 


sin 9 cos ( 0 — 0 



(Eq. 16) 


The exponential form is expanded in a Fourier series using the relationship 


exp j^iskr ^sinBcos ( 0— 0^)j = (i’^ (0 ~ 0 )) 

(Eq. 17) 


m= -“> 
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The point source Green's fimction then becomes 


4TTr 


00 

exp(iskr) ^ (skr^sine) j^exp ( im (0 - 0^ )) (Eq. 18) 


m=-<» 


As in the previous section, the net reaction forces on the acoustic media due 
to thrust and drag forces is represented as a distribution of dipoles over the 
rotor disk. The Green's functions representing the pressure contributions 
from thrust and drag dipole orientations are given as follows 


o z dr 


iskcos9exp (iskr) 
4TTr 


00 


E 

m= -00 


1 J 


m 


I skr sin9 ) exp 


(im (0-0^)) 
(Eq. 19) 


g0 = 


1 

r 90 


1 

r 


oo 

(exp (iskr) /4TTrj ^ ^ 


.m 

mi J 


m= -00 


m 


skr sin 
o 


9 1 exp(im( 9-9^ 
(Eq. 20) 


As before, the actual distribution of force over the rotor disk is modeled as 
acting about a ring of effective radius a. Using the previous expression for 
thrust distribution f^, the sound contribution from the thrust distribution is 
expressed as 





f^g r dr d0 
z o o o 


ikcos Q 
4TTr 


00 



s, m, q=-co 


2TT r 

/ / Pq 

o o 


J (skr sin9|exp 
m\ o / 


i (sB+q j 0^ 


— isuut + im 


(0-0^)]6(r-r^) 


r dr d0 
o o o 


(Eq. 21) 
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The 0Q integration yields zero terms except for m = sB + q, thereby eliminat- 
ing the summation over integer m,. The integration over the Dirac delta 
function finally provides the expression 



2TTaF° 

T 


cos 9 
4TTr 


^ ^^i(sB+q)+l 

S, 


a 6 J „ 

s q sB+q 


(ska sin 9) 


exp j^iskr+i (sB + qj 0-ismjt 


(Eq. 22) 


The derivation of an acoustic pressure component due to dipole drag forces 
proceeds in a similar fashion. The drag-force distribution is taken, as 
expressed in the following equation, to act at the same effective radius (a) 
as the thrust distribution 


CO 

s, q=-® ■ ' 


where F 


o Average Drag Force _ . D 

D Length at Circumference 2rra 


Utilizing the properties of the Green’s function 


2TT r 




Jq Jo 0 D o o o 


(Eq. 24) 


Performing the indicated integration yields 


= 


D 4TT 


^ X) ( = B+q)a^6^J^g^^(ska sin e 


s, q=-® 


X exp l^iskr+i (sB+q)— isujtj 


(Eq. 25) 


1 1 



I I null ■ I iini 


Note that the harmonic drag coefficients 6q, like the thrust coefficients Bq» 
are in general complex values. Without much loss in generality, some 
simplification is achieved by taking the phases of the drag harmonics to be 
identical with those of the thrust harmonics. Mathematically, 

6 = |6 I exp (iy ) (Eq. 26a) 

q q' ^ q' 

and 

= |6 I exp (-iY^) (Eq. 26b) 

The expression for the total acoustic pressure at an arbitrary far-field 
reference point is given by the instantaneous addition of the thrust and drag- 
derived component pressures 


P = 


+ P 


D 


(Eq. 27) 


Combining the expressions for Pq- and Pp) results in 


P=4Tfr5Z [-skfJP.|cos9+F„i6 

s, q=~ 


(sB+q )' 


t l I V, • — I I 

q D q a 


a J „ (ska sin 9 
s sB+q ' 


X exp l^iskr + i ( sB+q ) (0 + ^ ) + iYq ~ isU)t 


(Eq. 28) 


By expanding the series and regrouping terms, it can be shown that with the 
use of the expression (Z) = J_^ (-Z), the series for P can be changed to 
the following trigonometric form 


P 




2a 


<JU 

£[- 


skF^lB^lcosB + F^l 


(sB-q ) ~ 


J ^ (ska sin 0 
sB-q \ 


X 



(sB+q)~ 

a 


sB+q 


( ska sin 6) 



(Eq. 29) 
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where 




sk (r-ct) + (sB-q) (^ + ^ ) + 

tt 

sk(r-ct) + ( sB+q ) (0 + — ) + Y 

Z q 


(Eq. 29 ) cent. 


As discussed in References 3 and 6, the characteristics of Bessel functions 
related to the relative magnitude of the argument and the order may be used 
to make a simplifying approximation to the above equation. It is convenient 
to first express the argument of the Bessel function as 


skasinQ = sBMsinG 


where M = effective rotational Mach number 



(Eq. 30) 


For typical helicopter operations, this rotational Mach number is less than 
unity. The absolute magnitude of the argument is then always less than the 
integer value sB. For these conditions, the magnitude of the JsB+q functions 
is much smaller than the magnitude of the JsB-q functions. For the calcula- 
tions made in this study, the small contributions due to terms with a factor of 
JsB+q were neglected. 


The asymptotic properties of the Bessel functions were also used to set a 
limit on the number of terms included in the summation over the blade -load 
harmonic number q. The sum criteria recommended in Reference 3 is to 
include only the terms that satisfy the criteria 


sB (l-M ) < q < sB ( 1 + M) 


(Eq. 31) 


This criteria is strictly true only for a relatively flat blade -loading spectrum. 
For the calculations in this study, the lower limit of this criteria was re- 
moved. 


The expression for the time-varying pressure is finally written 


P = 


CO CD 


2a 


s = l q=0 



(Eq. 32) 
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X 


J^B q (sBMsinS ) sin j^sk (r-ct) + 


(sB-q)(0 + f)+ Yq 


(£q. 32) 


To compare predicted values with measured values, it is necessary to trans- 
form this expression to indicate mean-square harmonic amplitudes. The 
mathematical process of taking the time-average value of the squared pressure 
yields the expression 
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(Eq. 33) 


To make use of this general expression, it is necessary to specify each of 
the individual blade-load phases Yq as well as the harmonic amplitudes. 


An alternative, discussed in detail in Reference 3, is to assume that the 
phase amplitudes vary randomly in time. The time-averaged harmonic 
pressure amplitudes then become independent of the azimuthal variable 0 
and is expressed as 
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(Eq. 34) 
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3. 0 CALCULATED AND MEASURED RESULTS 


The approach of determining blade loads from on-axis acoustic measurements 
was investigated experimentally using both actual helicopter data and rotor 
tower test data. For each test series, the sound levels at a position on the 
rotor axis were spectrally analyzed to obtain a set of sound harmonic ampli- 
tudes. These sound harmonics were used with Equation 13 of Section 2, 1 to 
calculate rotor blade-load harmonic amplitudes. A least- squares curve fit 
was applied to the calculated blade-load harmonics to obtain an empirical 
relationship between harmonic number and amplitude. To test the utility of 
the empirical blade-load harmonic equations, they were used with the sim- 
plified random phase, rotational-noise prediction model (Equation 34) to cal- 
culate off-axis sound levels. These calculated sound levels were then com- 
pared with measured off-axis sound levels. The calculated and measured 
results are provided in the following report sections. 


3 . 1 Helicopter Test Series 

The helicopter noise data were obtained from a NASA-Langley noise measure- 
ment program in which the NASA Lunar Lander structure was used to position 
microphones about a helicopter in hover. Figure 3 is an illustration of the 
measurement configuration, showing the relative position of the helicopter to 
the microphone locations. Data were available from noise measurements of 
a four-bladed Hughes OH-6 helicopter and a two-bladed Bell 205 helicopter' 
shown in Figures 4 and 5, respectively. For the helicopter noise measure- 
ments presented in this report, each helicopter was maintained in a hover at 
an altitude of 24.4 meters (80 feet). In view of the height of the helicopter 
and microphone positions above ground, it was anticipated that although 
ground reflection interference of sound waves can occur, the effect will be a 
small perturbation and was consequently neglected. 


A block diagram of the instrumentation used to obtain the acoustic data is 
provided in Figure 6, The measurements were initially recorded on magnetic 
tape and subsequently processed using the SD 301C spectrum analyzer in the 
0 to 500-Hz range with a 1. 5-Hz bandwidth. The rms amplitudes of the 
rotational sound harmonics were read directly from the narrow band acoustic 
spectra. The microphone and signal amplification instrumentation provided 
a uniform frequency response from 20 Hz to lOkHz. Below 20 Hz, however, 
the system gain dropped sharply due primarily to the low-frequency response 
of the General Radio sound-level meter. The measured acoustic spectra ob- 
tained for the Bell helicopter tests showed that the first sound harmonic oc- 
curred at approximately 10 Hz. The amplitude of the first sound harmonic 
for these cases was consequently adjusted by 8 dB to account for the frequency 
response of the measurement systena. 
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The sornid harmonic amplitudes measured at 42. 7 meters ( 140 feet) directly 
overhead of the helicopter were used in Equation 13 to calculate blade-load 
harmonic amplitudes. For these calculations, the chord-loading coefficients 
a , were taken to be unity, and the effective radius taken to act at the 85% 
radial station. Since the helicopter was in a stationary hover, the total up- 
ward thrust was taken to be equal to the gross helicopter weight. 


The calculated blade-load harmonics for the OH-6 helicopter are shown in 
Figure 7, For the Bell 205 measurements, reliable data sets were available 
for two helicopter weight conditions, termed "low fuel" and "full fuel". The 
calculated blade-load harmonics for these two cases are shown in Figures 8 
and 9, As seen in Figures 7 through 9, the calculated blade-load harmonics 
tend to cluster about a straight line when plotted in a full log format. The 
straight line indicated in each of these figures was determined from a least- 
squares curve fit of the data points, and in the equation form shown represents 
an inverse power law relationship. In terms of decibel notation, the indicated 
decrease of blade-load harmonic amplitude with increasing harmonic number 
corresponds to a fall-off rate of approximately 6 dB/octave for the OH-6 
helicopter test case and approximately 8 dB/octave for the two Bell 205 heli- 
copter test cases. 


As a test of the empirical blade-load harmonic equations, they were used as 
input data in a prediction model to calculate off-axis sound levels. The ra- 
domized phase prediction model (Equation 34) was selected since the pro- 
cedure used to determine the blade-load harmonic amplitudes does not pl-o- 
vide any information about the loading harmonic phase angles. In addition, as 
discussed in References 2 and 3, the phase angle appears to be important 
only for the lowest order harmonics. Based on helicopter performance data, 
the ratio of total drag to total thrust was calculated to be 0. 104 for the OH-6, 
0. 096 for the low-fuel Bell 205, and 0, 106 for the full-fuel Bell 205 helicopter 
test cases. Therefore, in the prediction model, each blade-drag harmonic 
amplitude F^ 6 was taken to be 0, 1 of the corresponding blade-thrust har- 
monic amplitude F P . 

p T q 


The calculated directivity patterns for a sampling of the OH-6 rotational sound 
harmonic levels at a 52. 1-meter (171-foot) radius is shown in Figure 10. Also 
plotted in this figure are the measured sound harmonic levels normalized to a 
52. 1-meter (171-foot) distance by accounting for spherical spreading. The 
deviation between measured and calculated sound levels for the overhead 
microphone position is due to the errors in approximating the determined 
blade loads by the empirical power law equation. Figure 10 shows that the 
theoretical model for the OH-6 test configuration predicts a sharp drop in the 
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harmonic sound levels for elevation angles between TO'^and 90°. As shown, 
the measured harmonic levels at the 73°elevation are consistently higher than 
the calculated levels. The measured levels at the 36°elevation do, however, 
show scatter about the predicted values. These characteristics are also seen 
in the data for the "full-fuel" Bell 205 helicopter as shown in Figure 11. 


A direct comparison of calculated and measured sound harmonic amplitudes 
for a range of harmonic numbers is presented in Figures 12 through 14 for 
the OH-6, "full-fuel" Bell 205, and "low-fuel" Bell 205 test cases, respec- 
tively. The missing data points in these figures were cases in which a tail 
rotor sound harmonic partially masked a main rotor sound harmonic. 


3. 2 Rotor Tower Test Series 

In addition to the helicopter noise measurements, acoustic data were also 
obtained from a series of tests using the NASA- Langley Research Center 
helicopter rotor test facility pictured in Figure 15. The test series, identi- 
fied as rotor tower test 12-A, was conducted with a two-bladed type UH-ID 
rotor. A sketch of the rotor test facility showing the relative microphone 
locations is provided in Figure 16. The on-axis microphone was suspended 
directly above the rotor hub from a horizontal overhead cable. For the test 
series, the UH-ID rotor was operated at various pitch angles and rotational 
speeds. Data from six different rotor conditions were chosen for detailed 
study for this investigation. Using the measured on-axis sound harmonic 
amplitudes as data input, Equation 13 was used to calculate the blade-load 
harmonic values. The calculated values for the six data sets are shown in 
Figures 17 through 22. As shown on each of these figures, a least- s quares 
curve fit was performed to establish an empirical power law relationship for 
the blade-load harmonic values. These empirical relationships were then 
used with the rotational noise mathematical model to predict sound levels to 
be compared with the off-axis acoustic measurements. 


To provide an indication of the theoretical directivity pattern, calculations 
were made of a sampling of sound harmonics as a function of elevation angle 
for one of the test runs. These harmonic levels were calculated for a 76. 2- 
meter (250-foot) radius and are compared in Figure 23 with measured values 
adjusted for differences in spherical spreading. 


From the location of the off-axis microphones at 0 = 90° and 6= 99. 6°, the 
measured results were indicative of sound radiated in and near the plane of 
the rotor disk. Examination of the prediction model (Equation 34) shows 
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that for observation angles ©near 90°, the soiond levels will be controlled by 
the blade-drag harmonic amplitudes. For the sound harmonic calculations, 
as done for the helicopter cases, the blade-drag harmonic amplitudes were 
taken to be 0. 1 of the thrust harmonic amplitudes. Since the initial calcula- 
tions showed a general tendency to underestimate the measured sound levels, 
the calculation procedure was repeated with the harmonic drag-to-thrust 
ratio increased to 0. 2. A comparison of the calculated and measured sound 
harmonics for these cases is provided in Figures 24 to 29. For comparative 
purposes, the experimentally derived net drag-to-thrust ratio was also de- 
termined. Using the thrust and torque values measured on the rotor shaft, 
estimates of the net drag-to-thrust ratio were made by assuming the net drag 
force to act at the 85% radial station. For the rotor test facility cases con- 
sidered, these net drag-to-thrust ratios were in the range 0. 0625 to 0. 0645. 
Use of these values would result in a decrease of predicted sound harmonic 
levels of up to 4 dB. 


In view of the proximity of the rotor and microphone locations to the ground, 
it is likely that the ground reflected sound has had an effect on the measured 
data. A precise accounting of the effects of ground reflection would require 
details of the actual spatial distribution and acoustic directivity of the rota- 
tional noise source, as well as information about the acoustic impedance of 
the ground surface. Some indication of the effect of ground reflection can 
be seen, however, by comparing trends in the measured data with results 
calculated for the idealized configuration of an acoustic point source situated 
over an infinite rigid plane. A calculation of the interference pattern caused 
by ground reflections for this ideal case was made using the analytical pro- 
cedure detailed in Reference 7. 


The calculated interference patterns are shown in Figure 30 for sound har- 
monics based on a rotor rpm of 324. As seen in the top part of this figure, 
the in-plane microphone position is likely to exhibit the greatest interference 
effect. By comparison with Figures 26 through 29 for tests with similar 
rotor rpm, the measured in-plane data clearly shows the superimposed inter- 
ference characteristic, particularly near sound harmonic numbers of 3 and 10. 
For the microphone position 0. 8 meter (2. 5 feet) above ground, the first cal- 
culated interference dip occurs well above the 20th sound harmonic number. 
Figure 30 shows the effect of a rigid reflecting surface for this configuration 
would be to add nearly 6 dB to the direct radiated sound. For the case of the 
overhead microphone position, as inidcated at the bottom of Figure 30, the 
reflected sound causes and interference pattern with a variation of within 
+ 1 dB. 
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Because of the groimd cover, the interference effects caused by the actual 
ground reflection conditions at the helicopter rotor test facility are not likely 
to be as predonainate as the ideal case considered here. Therefore, no 
corrections for ground reflections were applied to the predicted data so as 
not to obscure the prediction trends by an oversimplified set of factors. 
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4.0 CONCLUSIONS 


The current study has provided a theoretical basis by which rotor-blade har- 
monic loading can be inferred by near-field as well as far-field acoustic 
measurements along the rotor axis. This analytic methodology has been 
applied to a variety of experimental cases involving noise measurements 
from two helicopter types as well as measurements taken during a series of 
static test facility runs. 


The derived blade-load harmonics for the helicopter and rotor test facility 
cases considered have shown a reasonable fit to a simple harmonic number 
power law expression. When expressed in decibel form, these power law ex- 
pressions provide decay rates of blade-load harmonic amplitudes with har- 
monic number in terms of dB/octave. A comparison of the blade-load har- 
monic decay rates determined from data in this report and decay rates dis- 
cussed by other investigators is shown in Table I. It was anticipated that air 
flow recirculation effects during operation of the rotor test facility might pos- 
sibly result in a harmonic blade load characteristic significantly different 
from that during a helicopter in hover. From a comparison of the experi- 
mental results summarized in Table I, it is apparent that the blade-load 
harmonic decay rates from both the helicopter and rotor test facility data fell 
within the range of 6 to 9 dB/octave. This is comparable to the 7 to 10 dB/oc- 
tave decay reported by Hosier and Ramakrishnan (Reference 2). As shown in 
Table I, the derived decay rates fall within the transition range between im- 
pulsive and nonimpulsive noise conditions reported by Bausch, Munch and 
Schlegel (Reference 8). Similarly, the derived results also fall between the 
values suggestive of blade slap and smooth inflow predicted by Lowson and 
Ollerhead (Reference 3). 


The acoustic predictions based on the empirically-derived blade loading power 
law relationships in conjiinction with the random-phase prediction model were 
of only limited success. For the helicopter cases considered, the predicted 
values were within about 5 decibels of the measured values for the measure- 
ment location at 36® from the vertical axis. For the measurement location at 
78®, the predicted sound harmonic levels were on the order of 10 dB lower 
than the measured values. This tendency of ander estimating the sound levels 
for positions near the plane of the rotor was also evident from the rotor test 
facility data. The predicted soxond harmonic for the 6= 90® and 0= 99. 6'- .ni- 
crophone positions were generally well below the measured values. 


There are several potential reasons for the lack of correlation between the 
measured and predicted data. From an experimental sti. i-.'dpoint, possible 
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explanations include acoustic measurement anomalies due to effects of ground 
reflection and the influence of micrometeorological conditions along the 
sound propagation path. Potential analytical explanations include prediction 
error due to the mathematical approximations inherent in the model and to 
the simplifications achieved by treating the loading forces in terms of a con- 
centrated ring source. An explanation which does correspond to the general 
data trends is that, in addition to the anticipated dipole rotational noise radi- 
ation, there was a significant contribution to the measured data from the 
monopole rotational noise source referred to as blade-thickness noise. As 
discussed in References 9 and 10, blade-thickness noise is caused by the peri- 
odic displacement of fluid due to passage of the rotor blade. Although blade- 
thickness noise has generally been recognized as a source of less importance 
for subsonic rotor tip speeds, the sound directivity pattern does peak in the 
plane of the rotor disk. 
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TABLE I 


Blade-Load Harmonic 


Data Source Fall-off Rate dB /Octave 

Remarks 

Figure 7 

6 

OH -6 in hover 

Figures 8 and 9 

8 

Bell 205 in hover 

Figures 17 to 22 

6. 5 to 9 

Rotor test facility with 
UH 1-D blade 

Hosier and Ramakrishnan 
(Reference 2 ) 

7 to 10 

Rotor test facility with 
four-bladed hingeless 
rotor 

Bausch, Munch and Schlegel 
(Reference 8 ) 

8. 5 to 9. 5 

Nonimpulsive noise 
conditions 

Bausch, Munch and Schlegel 
(Reference 8 ) 

3. 5 

Impulsive noise 
conditions 

Lowson and Ollerhead 
(Reference 3 ) 

15 

Predicted for smooth 
inflow 

Lowson and Ollerhead 
(Reference 3 ) 

0 to 6 

Suggested for blade 
slap conditions 
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Figure 6. Instrumentation Block Diagram. 
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Figure 7. Calculated Blade-Load Harmonics; 
Lunar Lander Test Series, 
Hughes OH- 6 
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Figures. Calculated Blade-Load Harmonics; 
Lunar Lander Test Series, 

Bell 205 (Low Fuel) 
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Figure 10. Calculated Directivity Pattern of Sound 

Harmonics with Comparison to Measured 
Values for Hughes OH-6, Lunar Lander Series 
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Figure 11 . Calculated Directivity Pattern of So\ond 

Harmonics with Comparison to Measured 

Values for Bell 205i(Full Fuel), Lunar Lander Series 
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Figure 14. Comparison of Measured and Calculated 

Sound Harmonics for Bell 205, (Full Fuel); 
Lunar Lander Test Series 
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Figure 15. NASA Langley Research Center 
Helicopter Rotor Test Facility 
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Figure 16. Microphone Locations for Rotor Tower Test 12-A 
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Figure 23., Calculated Directivity Pattern of Sound 

Harmonics with Comparison to Measured 
Values for Rotor Tower Test 12-A, Run 112 


47 







Sound Pressure Level, dB Sound Pressure Level, dB 



,48 


Figure 24., Comparison of Measured and Calculated Sound 
Harmonics for Rotor Tower Test 12-A, Run 104 




Sound Pressure Level, dB Sound Pressure Level, dB 





Sound Pressure Level, dB Sound Pressure Level, dB 


100 , 


In Plane Microphone 
{ 0 =90°) 


• Measured Values 
- Calculated Values 



5ol__l I l_l I I I I I l_l I I I Lj I I L 



Figure 26 . 1 Comparison of Measured and Calculated Sound 
Harmonics for Rotor Tower Test 12-A, Run 112 
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Figure 27. Comparison of Measured and Calculated Sound 
Harmonics for Rotor Tower Test 12-A, Run 114 


51 



Sound Pressure Level, dB 


• Measured Values 
- Calculated Values 



Figure 28. Comparison of Measured and Calculated Sound 
Harmonics for Rotor Tower Test 12-A, Run 122 
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Figure 29. Comparison of Measured and Calculated Sound 
Harmonics for Rotor Tower Test 12-A, Run 126 
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